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Finite T electroweak phase transition on the lattice^ 

K. Rummukainen a * 

a Indiana University, Department of Physics, Swain Hall- West 117, Bloomington, IN 47405, USA 

This talk reviews recent lattice results on the high T electroweak phase transition. A remarkably accurate 
picture emerges: a) the transition is of first order for m_ff<80GeV and vanishes for larger m,H\ b) transition 
temperature, latent heat and interface tension are known, as well as c) the properties of the broken and symmetric 
phases. New developments in the sphaleron rate calculations are discussed. 



1. INTRODUCTION 

At high temperatures the spontaneously bro- 
ken symmetry of the electroweak (EW) sector 
of the standard model is restored Q. How this 
restoration occurs is crucial for understanding the 
effects of EW physics to the baryon number of 
the Universe, and for the viability of EW baryo- 
genesis ||. The essential questions are: what 
is the order of the transition, or is there only a 
smooth cross-over? What are the values of the 
transition temperature T c , latent heat L, inter- 
face tension a and the discontinuity of the Higgs 
condensate? Very important are also the prop- 
erties of the phases close to T c : the equation 
of state, screening lengths, and phase metasta- 
bility ranges. For baryogenesis scenarios, a cru- 
cial quantity is the sphaleron rate, i.e. the rate 
of the baryon number fluctuation. These quanti- 
ties are parametrized by the still unknown Higgs 
mass ran (experimentally m#>64GeV). All of 
the questions above have been addressed with 
lattice simulations; and happily, more often than 
not, definite answers have been found. 

Given the success of the EW perturbation the- 
ory at T — 0, it is natural to apply it at finite T: 
indeed, the effective potential has been calculated 
up to 2-loop level H|). However, it has become 
clear that the perturbation theory in gauge the- 
ories at high T is intrinsically unreliable due to 
infrared divergences |J. At low T the perturba- 
tive expansion is regulated by the large value of 
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the Higgs condensate v: the expansion parameter 
is ~ g 2 T/[irmw)- When T increases, v becomes 
smaller, and in the symmetric phase the expan- 
sion cannot be controlled even with resummation 
techniques. Clearly, a non-perturbative approach 
is required. 

In comparison with the QCD phase transition, 
until recently the finite T EW transition had 
not been extensively studied. Since the non- 
perturbative effects are expected to be mainly due 
to the SU(2) gauge fields, the studies have con- 
centrated on the SU(2) gauge + Higgs model. A 
direct way to study the finite T physics is to per- 
form 4D Euclidean lattice simulations; the stan- 
dard formalism was set up almost a decade ago 
HQ, but the first results close to the physical 
weak coupling parameter values were published 
only in 1992 g. The development of the di- 
mensionally reduced 3 dimensional effective for- 
malism, initiated in || and completed in p|,|l0||, 
was an important milestone theoretically and es- 
pecially for practical simulations. Until Lattice 
'95, Higgs mass ranges vhh = 18-49 GeV had 
been investigated with 4D simulations 111 
and 35-80 GeV with the 3D formalism g|4 
(for earlier reviews, see ^7|-|l9|). In this confer- 
ence new 4D results with run 5^102 GeV were pre- 
sented by the DESY group @J§ and Y. Aoki 
fH, and 3D results with m H < 180 GeV by Ka- 
jantie et al.|J,||], Karsch et aU|lil> Giirtler 
et al.[^8| and Philipsen et al.p9[! New studies 
of the sphaleron rate were reported by Tang and 

smit iqi. 
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2. THE SU(2)-HIGGS MODEL 

Since the essential physics of the EW phase 
transition is expected to arise from SU(2) gauge 
fields and the Higgs field, let us for the moment 
neglect SU(3) and U(l) components of the gauge 
fields and the fermions entirely. The SU(2)-Higgs 
Euclidean Lagrangian is 

l = \f; u f; u + (d^)\d^) 

- m 2 <j>U + H^<f>) 2 , (1) 

where (j> = (0 + ,0°) is the Higgs doublet. In the 
following we discuss the derivation of the effective 
3D action from ([I]). 

2.1. 3D effective action: why and how 

Because of the relatively small value of the cou- 
pling g, the EW theory at high temperatures has 
a very wide range of mass scales (~ inverse screen- 
ing lengths): 

T»m I? « 5 T» g 2 T,m w (T),m H {T). (2) 

For example, the simulations have shown that if 
m H w 60GeV, then m H (T c ) < 0.1T C . On an 
Euclidean finite T system the imaginary time is 
restricted to the interval < r < 1/T; thus, 
when this system is latticized, in order to avoid 
finite size effects the ratio of the spatial and tem- 
poral extents of the lattice should be at least 
L S /L T > 30-40. This is a very punishing require- 
ment. 

The extreme 'flatness' of the geometry already 
suggests that the essential features of the sys- 
tem might be described by an effective 3D the- 
ory. Because of the periodic boundary condi- 
tions, when the bosonic fields are expanded in 
Fourier modes the inverse propagator becomes 
[fc 2 + m§ + (27mT) 2 ]. If T is large compared to 
the other relevant mass scales of the system, the 
non-static Matsubara modes with n = ±1, ±2, . . . 
acquire a heavy mass lixnT and can be integrated 
over. What remains is an effective 3D theory of 
the static (n = 0) modes. If the original 4D the- 
ory has fermions, then, because of the antiperi- 
odicity in t, all fermionic modes become mas- 
sive with m — (2n + l)nT and can be integrated 
over. The effect of the fcrmion fields is only to 
adjust the parameters of the effective 3D bosonic 



theory. The derivation of the 3D effective the- 
ory, dimensional reduction (DR), introduced in 
| 3l|. was fully developed for the EW theory in 

10 it- 

The integration over the non-static modes can 
be performed perturbatively. This is possible if 

• 9 2 , A < I 

• T> mg(T c ), the relevant mass scales at T c . 
For the EW transition, these conditions are met 

when 30<m# <240 GeV: the lower bound comes 
from the requirement that m#(T), mw{T) must 
be <C T in the vicinity of T c , and the upper bound 
from the fact that for large m# the EW theory 
becomes strongly coupled. Note that there are no 
IR problems in the derivation of the effective ac- 
tion: the IR modes are inherently 3-dimensional 
and are not integrated over. Indeed, the effec- 
tive theory retains all IR divergences of the orig- 
inal theory! Moreover, the perturbative DR does 
not require that the finite T perturbation the- 
ory should be applicable in general: the crite- 
rion for finite T perturbation theory to work is 
g 2 T/niQ <C 1, which is not satisfied at high T. 

Starting from eq. ([!]) the effective theory can 
be defined by the action Ss[Aq(x), Af(x), 4>(x)], 
a 3D SU(2) gauge + adjoint Higgs + fundamen- 
tal Higgs theory, where the coefficients of the ac- 
tion depend on the original 4D action coefficients. 
The adjoint Higgs Aq is the remnant of the time- 
like component of the 4D gauge field. This ac- 
tion was used in [p|JT4|] to simulate mjj = 35 and 
80 GeV systems. Further, for the EW theory, one 
observes that the Debye mass mjj = \Jb/d>gT is 
large, and the field A can also be integrated over. 
The resulting action 

Ss[A h( t>\ = J ^x[\F^F^ + {Di4^iPi<i>) 

+ m^ + A 3 (<^) 2 ] (3) 

has been the 'workhorse' in all recent 3D sim- 
ulations |p4] , p8| , p6| . It is similar in form to the 
original 4D action (|l|); however, now the fields 
and couplings have dimensions [0] = GeV 1 / 2 , 
[g 2 ] = [A3] = GeV. In Q the action was fur- 
ther simplified to an 0(4)-symmetric scalar the- 
ory by integrating over the gauge fields. However, 
this action failed to describe the phase transition 
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correctly, indicating the essential role of the mag- 
netic sector of the gauge fields. 

The theory is now uniquely fixed by three pa- 
rameters, the gauge coupling g\ and two dimen- 
sionless numbers 



2\2 



x = h/g 3 , V = m 3 (g 3 )/(g 3 ) 



(4) 



The action is superrenormalizable: the cou- 
plings g 2 and A3 do not run (in MS) , and m§ has 
only linear and log-divergences at 1- and 2-loop 
levels. This property makes the relation between 
lattice action and continuum parameters (0) par- 



ticularly transparent, as discussed in section 2.3 
The action of form (^) can be derived already 
with 1-loop DR. However, to fully utilize the su- 
perrenormalizability better accuracy is required, 
and in ||,[Io]j24] DR is performed with Green's 
junction matching : one starts from a general 
3D superrenormalizable action and matches all 
2- and 4-point Green's functions to the static 4D 
Green's functions of the original theory. To can- 
cel large logarithms (in MS), the 4D couplings are 
run to scale /i = 47rTe~ 7 « IT by the standard 
4D j3 functions. The matching is done to a consis- 
tent accuracy in g 2 , A. The action (^) can provide 
the relative accuracy 5G/G ~ 0(g 3 ). To go be- 
yond this would require the inclusion of 6-dim. 
operators while giving up superrenormalizability. 
By calculating the effects of the excluded higher 
dimensional operators the accuracy of the effec- 
tive action can be estimated |lC| ]. By construc- 
tion, the Green's function matching avoids the 
non-localities inherent in the standard integra- 
tion and 'block transformation' -type approaches 
to effective actions Q . 

The DR process provides us with the essen- 
tial relations between the 3D and 4D parame- 
ters. There is a large class of 4D theories which 
map into the 3D SU(2)+Higgs theory. Since the 
fermions do not 'survive' DR, this class includes 
the 4D theory of SU(2) + Higgs + fermions and 
the minimal standard model (MSM), where the 
effect of U(l) gauge field can be estimated per- 
turbatively. The mapping [g 2 1 x,y) <->4D param- 
eters for these theories has been worked out in 
detail in |l(|. Conversely, a single 3D simulation 
can yield physical results for the whole class of 
4D theories. Recently DR has been worked out 



for MSSM by several authors |]35j. 

For the 4D SU(2)-Higgs theory the relation 
4D^3D is MM 



9l = 



0.44015T* (5) 

x = -0.00550 + 0.12622/i 2 (6) 
y = 0.39818 + 0.15545/i 2 

- 0.00190/i 4 - 2.58088(m^/T*) 2 , (7) 

where h = m* H /mw, raw — 80.6 GeV and 



A = \9 2 { m *Hl m w)' 2 



5=3' A = g9"[m' H /mw)" ( 8 ) 

(The authors of |^,|8| use somewhat different 
conventions.) The notation m* Hl T* is used to re- 
mind that these are not the physical T or the pole 
itih', for the SU(2)-Higgs theory without fermions 
the difference is small. These parameters are used 
extensively to present the results of 3D simula- 
tions. 

2.2. Lattice action in 4D 

The lattice action is conventionally written as 

(3 G J2 

fe^^$tu i; ^ + .) (9) 

x; fi 

£ ^Tr < $ x + 0R £ [ Jtt < $ x 1] 2 



S = 



where $ = RxV, V £ SU(2), and R 2 = ±Tr$t$. 
The Higgs field $ has SU(2) gaugc ® SU(2) global 
symmetry. 

The essential question is now the relation of the 
lattice parameters to continuum physics. In the 
4D formalism this is done non-perturbatively, re- 
lying only on the measurements of physical quan- 
tities. This comes at a cost: simulations at T = 
are needed in order to set the scale. At tree level, 
the relation is 



0G = 

(ma) 2 
A = 



4/3 2 

(1 - 4/3 ff - 20 r )/0h 



(10) 

(11) 
(12) 



Strictly speaking, for the 4D SU(2)-Higgs the 
a = limit cannot be reached because of triv- 
iality. The term "continuum limit" in this case 



4 



means reaching a region where the cut-off ef- 
fects become negligible: physics remains invari- 
ant when moving along the renormalization group 
trajectories, or constant physics curves (CPC). 

The most detailed scaling study so far has been 
performed by the DESY group |2(| using lattices 
with N T = \/(aT) = 2-5. A practical way to 
determine CPCs is to use the transition itself as 
follows @,|0): 

(a) To be close to the desired physical point 
g 2 ~ 0.5, m H w 34GeV, couplings (3 G = 8, 
Pr = 0.0003 are chosen for N T = 2 simulations. 
mw is fixed to 80 GeV. The value of Ph is tuned 
until the transition coupling Ph, c is found. 

(b) Using these couplings one performs T = 
(N T >N S ) simulations. Higgs and W masses and 
the static potential (Wilson loops) are measured; 
from the static potential the renormalized gauge 
coupling g 2 R can be extracted. This also gives 
-Vh = 9R/^{jn H /m w ) 2 . 

(c) Using the known continuum 1-loop RG- 
equations for Pq and /3r, one runs along CPC 
from N T -> N T + i ( a _> a N T /(N T + 1)). The 
step (a) is then repeated with the new couplings. 

Good scaling now means that dimensionless 
physical quantities g R , \r, T c /mw remain in- 
variant. This is surprisingly well satisfied already 
when N T = 2 — > 3, in strong contrast to QCD or 
even pure gauge SU(N) phase transitions, where 
scaling violations are still seen at N T = 6. Es- 
pecially the T = quantities g R and mn/mw 
do not show systematic a-dependence beyond the 
statistical errors. In this case g 2 R w 0.585 and 
mn/inw ~ 0.422, close to the tree level value. 
The good scaling behaviour indicates that the 
modes constant in r are the dominant ones. 

In the physically allowed range tuh ^64 GeV, 
the inequality T ^> to#(T), mw(T) makes the 
4D formalism very unwieldy. An appealing rem- 
edy for this is to use asymmetric lattice spac- 
ings a s = £a T . The kinetic part of eq. (|^) splits 
into temporal and spatial parts with couplings 
(3q, (3q and P T H , P S H . The couplings are related 
to the asymmetry through relations Pl/Pf = 
These have been evaluated perturbatively 
to order 0(g 2 , A) by requiring isotropy in the W 
and Higgs correlations P^| . In non-perturbative 
tests these relations were reproduced very accu- 



rately [g2|. 

2.3. Lattice action in 3D 

The 3D SU(2)-fundamental Higgs lattice action 
is similar in form to the 4D action in eq. ([}]), ex- 
cept now the indices are limited to values 1-3. 
The essential difference between the 4D and 3D 
formalisms becomes evident when we look at the 
derivation of the CPCs. The crucial point is the 
superrenormalizability of the 3D action (Q): g\ 
and A3 do not run, and rn^(g^) has only 1- and 
2-loop divergences. Calculating the relevant dia- 
grams both in the continuum and on the lattice 
one obtains the relations |[l(XB4l 



4 

4 Wg = ~W 



1 

fa 



2xfa\ , 3S/3 G 



16tt 2 



51 
16 



9x- 12x 2 



(1 



32tt 



(13) 
(14) 

h4a;) 



5.0 + 5.2a; 



(15) 



Eq. ([l5]) depends on several constants arising from 
lattice perturbation theory: £ = 3.17591, ( = 
0.09, and the two numbers 5.0 and 5.2, specific 
for SU(2)-Higgs theory and calculated in p3|. 

For fixed continuum parameters (g^x, y) 
eqs. (|l3|-|l5|) define the CPC in the space of 
(Pg, fa, fa) when lattice spacing a is varied. 
Also note that in 3D Pa oc 1/a. In contrast to 
the 4D case the continuum limit is well defined in 
3D. The equations above have relative accuracy 
0{a~ l ), so that the relation continuum-f->lattice 
becomes exact when a — > 0. However, in practice 
the finite a effects have been observed to be small. 

The Bielefeld group takes a different 

philosophy to 3D effective theories: they do not 
utilize the superrenormalizability of the 3D ac- 
tion, but consider that the most natural approach 
is to fix the cutoff scale to be of the same order 
of magnitude than the physical scales. While this 
introduces differences oc a, the smallness of the 
finite a effects makes most of the results compa- 
rable to the a — ► results. 
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3. PHASE TRANSITION 

3.1. Phase diagram 

In previous Lattice meetings 

HQ results 

with ma from 18 to 80 GeV were reported. The 
transition was seen to be quite strongly 1st order 
at small m# , and to weaken rapidly with increas- 
ing mji ■ When m# = 80 GeV the resolution was 
not good enough to distinguish the order of the 
transition. Since this mass region is expected to 
be physically relevant, it is essential to clarify the 
situation at larger ran. 

An important point in understanding the phase 
structure is the observation that the EW theory 
does not have a true gauge invariant order pa- 
rameter which would distinguish the symmetric 
and Higgs phases. Indeed, it has been shown an- 
alytically that in the SU(2)-Higgs lattice system 
where the Higgs length is fixed the Higgs and the 
confined phases are analytically connected ^] in 
the strong coupling regime; this was also observed 
in the early simulations ■ 

A study of one-loop Schwinger-Dyson equa- 
tions 1 38 1 argues in favour of the end of the 
1st order transitions at m# ~ 100 GeV, after 
which only a smooth cross-over remains. This is 
certainly consistent with the observations above. 
However, the result relies on the applicability of 
the perturbation theory which is known to break 
down at m# ~ 80 GeV. On the other hand, e- 
expansion techniques |59| predict that a weak 1st 
order transition remains even with large run . Due 
to the lack of an order parameter, it is not likely 
that the 1st order transition turns into a line of 
2nd order transitions. 

The situation has been addressed by the recent 
3D and 4D Q simulations, which indi- 

cate that the transition turns into a smooth cross- 
over at ms ~ 80 GeV. Kajantie et al.|^5j utilized 
the finite size behaviour of the ^<j> susceptibility 
(in 3D notation): 

X = glVMU - (<^)) 2 ) (16) 

where 1 9\ = /3 g I3h(R 2 )/8 + const, and 

V(g^) 3 = (4A//3 G ) 3 . For each V the tempera- 
ture is adjusted until the maximum value of \ 
is found. According to the finite size scaling, 
in 1st order transitions Xmax °c V, in 2nd or- 
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Figure 1. X (T) for m* H = 60 and 120 GeV around 
the maximum |25| (Note the different y-axes!). 



der transitions x max oc V 1 ^ 3 , where 7 is a criti- 
cal exponent, and if there is no transition x max 
approaches a constant value. In Fig. |l| x(T) is 
shown for m* H = 60 and 120 GeV. The difference 
is striking: the quantities Xmax/V (60 GeV) and 
X (120 GeV) approach constant values, consistent 
with a 1st order transition and no transition, re- 
spectively. Xmax for m* H = 35-180 GeV are plot- 
ted against V in Fig.||. Scaling is well satisfied: 
points with different (3g fall on the same curve. 
To emphasize the approach to the asymptotic V 1 , 
V -lines, a simple mean field model has been fit- 
ted to the data; the results of the fits are shown 
as continuous curves. 

Similar behaviour was reported by Y. Aoki P3[ , 
using 4D rriH — 47-102 GeV simulations, and by 
the Bielefeld group §7j|, using m H = 60-100 GeV 
in 3D. Interestingly, also in finite T U(l)-Higgs 
theory an endpoint of the 1st order phase transi- 
tions has been observed fijfl . 

A more detailed study of the properties of 
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Figure 2. Xmax for different m* H as a function of 
V. The continuous lines are mean field fits, the 
dashed line corresponds to the mean field critical 
exponent. 



the critical point was performed in p7|| . In the 
mjj = 80GeV case i71h(T) was measured around 
the transition in the symmetric and broken phases 
and fitted to the ansatz 




m H oc \(3 G - /3 G ,c\ l 



(17) 



(here (3g is adjusted while keeping (3h, (3r con- 
stant), with the result v = 0.49(2) and 0.31(1) in 
the symmetric and broken phases, respectively. 
The different indices could indicate a tricriti- 
cal nature for the endpoint; on the other hand, 
80GeV is likely not the exact value of itih,c- In- 
deed, utilizing the analysis of the Lee- Yang zeroes 
for rriH = 80 and 100 GeV it was estimated that 
mH,c ~ 77 GeV. 

3.2. T c and metastability 

To accurately locate T c both the finite size and 
finite lattice spacing effects have to be addressed. 
In 3D simulations, for fixed (3g = 4/(g|a), x and 
V, one adjusts y in order to find the pseudo- 
critical point (Ph and (3r are given by eqs. (jlj- 
|l5|)). Through eqs. (||-|7|) this corresponds exactly 
to adjusting T* while keeping m* H constant. 
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Figure 3. Top: the V — > oo extrapolation of the 
pseudocritical coupling 0h,c- Bottom: the contin- 
uum limit (1/Pg = 0) extrapolation of T c (open 
symbols). Also shown here is the metastability 
range (filled symbols) pi ]. 



With finite V, there are several non-equivalent 
methods to locate the pseudocritical point: the 
maximum of the order parameter susceptibility, 
the minimum of the Binder cumulant and the 
"equal weight" and "equal height" points of the 
order parameter distributions |24| , ^8|] , In the 
V — > oo limit these all extrapolate to the same 
value, as shown in top part of Fig. ^. Multicanon- 
ical simulations and histogram reweighting are 
commonly used in the analysis. This is repeated 
for 2-4 values of (3c while x is kept constant. The 
continuum limit is obtained by extrapolating in 
1/Pg- This is shown in Fig. || for m* H — 60 GeV; 
in this case the range in /3,3-values (5-20) is so 
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Table 1 

T* and T c for physical SU(2)+Higgs and physical 
SU(2)+Higgs+fermions -theories. m* H = 35-70 
is from @, 72.18 from || (70GeV in authors' 
notation). All units are in GeV. 



m* H 35 60 


70 


72.18 


T* 92.64(7) 138.38(5) 


154.5(1) 


157.74(5) 


T* pert 93.3 140.3 


157.2 


160.9 


SU(2)+Higgs 






m H 29.1 54.4 


64.3 


66.5 


T c 76.8 132.6 


151.2 


154.7 


SU(2)+Higgs+fermions (m top = 175 GeV) 


iriH 51.2 


68.0 


69.4 


T c 89.8 


105.8 


107.2 



Table 2 

T c from 4D SU(2)-Higgs simulations. 



mn/GeV T c /mw 



16 

34 
48 



0.464(2) 

0.910(16) 

1.153(16) 



3 

oo 
3 



12] 
|2C] 
12] 



large that the subleading behaviour is seen. Nev- 
ertheless, it should be noted that the T* varies 
very little across the extrapolation: the curvature 
is seen only because of the very small statistical 
errors. 

The phase metastability range is also shown 
in Fig.[| This has been obtained as follows: at 
y c , the order parameter histograms have a 2-peak 
structure. The histograms are reweighted off y c 
until the "shoulder" of one of the peaks van- 
ishes; in terms of the constrained effective poten- 
tial this corresponds to the temperature at which 
the metastability of the phase vanishes |24j . 

The 3D effective theory describe s a whole class 
of 4D theories, as discussed in sec. |2~T| . In table @ 
the physical T c is shown for 4D SU(2)+Higgs and 
SU(2)+Higgs+fermions -theories, in addition to 
the "bare" T* -values. The fermion content is the 
same as in the MSM, with mt op = 175 GeV. Note 
that the physical ran is different from m* H , and 



that the 



•h 



— 35 GeV -set does not correspond 



to any physical fermion theory. 

Locating the critical couplings in 4D simula- 
tions is essential for determining CPC, as dis- 
cussed in sec. 2.2 . The methods above can be used 
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0.95 - 
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■ 3D 



v N=2 



20.0 30.0 40.0 50.0 60.0 70.0 80.0 
mJGe\ 

Figure 4. T c /TP° rt from 4D @ and 3D @ sim- 
ulations. The open symbols show the extrapo- 
lation to the continuum limit (only m* H = 60 for 
3D) and have been shifted horizontally for clarity. 



also in 4D to locate the critical couplings; in ad- 
dition the DESY group has used methods which 
rely on the coexistence of two phases in long 
cylindical lattices [ p"2"|]2"c| ] : in the "constrained" 
method the order parameter is restricted to a nar- 
row range between the pure phase values, enforc- 
ing the system to reside in a mixed phase. The 
coupling is then tuned so that the distribution 
in this region is horizontal. This is equivalent to 
the condition that the 'flat' part of the 2-peak 
histogram is horizontal, but requires much less 
cpu-time than the full histogram calculation. In 
the "2-coupling" method the system is split into 
2 halves, and the critical coupling is bracketed 
by tuning the couplings in the two subvolumes 
individually so that the 2-phase configuration is 
maintained. 

In table || the ratio T c /m\y is shown for dif- 
ferent m,H- For the 34 GeV case results from 
N T = 2-5 have been used to extrapolate a — > 
quadratically. 

A direct comparison between 3D and 4D simu- 
lations at this stage is not straightforward, due to 
very different connections to continuum physics. 
The 3D simulations have used g 2 (7T) rj 0.444, 
whereas in 4D simulations g 2 B w 0.58. The re- 
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lations {g 2 ,x,y) ^4D quantities allows one to 
adjust g 2 ; however, then also the physical run 
changes. Currently there are no simulations 
which would correspond to the same physical sit- 
uation (however, see {H)). Nevertheless, it is 
straightforward to compare the results individ- 
ually to 2-loop perturbative results. In Fig.|| 
T c /TP crt is plotted, the squares correspond to 3D 
m* H = 35, 60, 70 and 80GeV cases, circles to 4D 
m,H = 34 GeV. Also shown is the approach to the 
continuum limit; for 3D the points are the same 
as shown in Fig. ||. A striking feature is the small- 
ness of the statistical errors in the 3D simulations. 
This is due to the rigorous nature of CPCs in 3D: 
the errors in T c are directly translated from the 
statistical errors of critical couplings. In 4D the 
errors accumulate through T = mass measure- 
ments. For the bare values of the critical cou- 
plings the accuracy is comparable, up to 6-7 dec- 
imal places. The finite a effects are seen to be 
small in both cases, and T c /T^ crt is consistently 
smaller than 1. Given the difference in the values 
(and renormalization schemas!) of g 2 and the 4D 
statistical errors the results can not be considered 
inconsistent. 

3.3. Interface tension and latent heat 

The tension of the interface a between the 2 
coexisting phases and the latent heat L of the 
transition are primary quantities characterizing 
the strength of the transition. The now ubiqui- 
tous histogram method [j42]| is the most accurate 
but numerically very demanding way to measure 



a: for each lattice 3-volume V 



- „3r.2 



L L z , where 



L < L z is assumed, one measures p max and p nimi 
the probability distribution maximum and the 
minimum between the peaks at T c . a is now ob- 
tained at the V — ► oo limit of the expression 
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(18) 



where G = ln3 for L x = L z and for Lj « Lj, 
and c is a constant. In addition, a — > limit 
should be taken. Eq. ( Jlq) h as been used both in 
3D and 4D simulations |J,||f||||. The DESY 
group |l2]]l3|j2l| , p2| ] has also used the 2- coupling 
integration method ]43| ] and the tunnelling cor- 
relation length MM analysis. These methods are 
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Figure 5. The interface tension (top) and latent 
heat (bottom) from 4D and 3D simulations. 



not quite as demanding computationally as the 
histogram method. 

In Fig.|| a is shown from several calculations. 
Both the 3D calculations shown perform a — > 
extrapolation. Note the dramatic decrease in a 
when rriu increases. Considering the difficulties 
in measuring a reliably the agreement must be 
considered to be good; the only point which dis- 
agrees somewhat with the trend is the 80 GeV 
asymmetric lattice (a T = as/4:, N T = 2) point 
presented in this conference p2|| . Perturbatively, 
a can be calculated with any reliability only when 
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the transition is strong: indeed, at very small run 
the agreement between the lattice and perturba- 
tive results is fair, but already at ~ 60 GeV the 
lattice results are a factor of 4-5 smaller. 

The latent heat L is directly related to the jump 
of {4>'4>) at the transition through the Clausius- 
Clapeyron equation: 

L/T^^A^U). (19) 

In Fig.|| L is shown from 3D and 4D simulations; 
in 3D the extrapolation to the continuum limit 
has been performed. Both in 3D and 4D also 
alternative methods to eq. have been used for 
measuring L, with practically unchanged results. 
Within errors, the results are consistent with the 
2-loop perturbation theory. 

3.4. Interaction measure 

The interaction measure 5 = e — 3p character- 
izes the deviation of the system from the mass- 
less ideal gas behaviour. This has been measured 
by the DESY group in 4D at m fl = 34 GeV 
p0[ ; in 3D S, e and p are not readily accessi- 
ble. When 0.5 < T/T c < 1 in the broken phase, 
(5/T 4 w 0.6, and at T c it jumps discontinuously 
to w 0.9. When T increases further 8/T 4 falls 
rapidly, and at T a 2T C it is ~ 0. This qualita- 
tively agrees with the perturbation theory in the 
broken phase, where it is applicable. 

4. BROKEN AND SYMMETRIC PHASE 

Deep in the broken phase perturbative analysis 
is well controlled. Lattice studies at T<T C can 
reveal the accuracy and the eventual failure of 
the perturbation theory. Indeed, as mentioned 
in the introduction, the expansion parameter is 
~ 5l/( 7rm M / ) 5 where mw(v) formally when 
the symmetry is restored. 

In refs. @,^8| the behaviour of 4>) is studied 
with 3D simulations. The general agreement be- 
tween the lattice and continuum 2-loop perturba- 
tive results is good, progressively becoming worse 
when T c is approached. However, since (4>^(jj) 
can be determined very accurately in the broken 
phase (relative error ~ 10~ 3 ), deviations are seen 
even relatively deep in the broken phase. This 
information was used in [E4fl to infer the value of 




Figure 6. mn and mw in the broken phase at 
m* H = 70 GeV compared to the perturbation the- 
ory (Giirtler et al.||). 

the so far uncomputed 3-loop term in the effec- 
tive potential, and to verify that it is linear in cf>. 
Even with this correction the effective potential 
fails at the transition point. 

The Higgs and W screening masses are mea- 
sured with the scalar (0 ++ ) and vector (1~~) op- 
erators of type 



S x = Tr[$t$ 3 



V a - 



Tr [T a & x U x ,& x+ ~ % 



(20) 
(21) 



where r° is a Pauli matrix. The masses (inverse 
correlation lengths) are extracted from the cor- 
relation functions (S x S y ) and (V x V y ); in the V- 
correlations only the diagonal part survives. To 
improve the signal the operators can be smeared 
or blocked in various ways. 

In Fig. | the measured values of mg and mw 
are compared against 1-loop perturbative values 
(with and without the wave function renormal- 
ization) for m* H = 70 GeV |28|. Deep in the bro- 
ken phase the results agree within errors; how- 
ever, closer to the transition deviations appear. 
Note that perturbatively the transition occurs 
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0, whereas in the actual simulations 
Similar behaviour has been observed 
15 



in other simulations |24 2 

The symmetric phase vector and scalar masses 
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Figure 7. Vector (filled) and scalar (open) masses 
in the symmetric phase |2q] . 



measured with operators fl20f]21[ ), are shown in 
Fig.0 @. Both masses increase when T (m|) 
increases. An interesting observation is that the 
m* H = 35 and 70 GeV masses are equal - this is 
reasonable, because in the symmetric phase the 
\(<jJ (j)) 2 -term in the action (|^) is very small. 
Comparing Figs.^j and ^ we note that there is 
a small but clearly discernible discontinuity of 
the masses at the transition: in the symmet- 
ric phase the masses become larger. Again, 
similar behaviour has been observed by other 
groups fHHJ. 

Philipsen et al.[^9| used various levels of block- 
ing of the operators (|2(i|-f2"l|) , and measured the 
full correlation matrix between different block- 
ing levels. By performing an eigenstate analysis 
the ground state and a few lowest exited states 
could be distinguished. The ground states in the 
vector and scalar channels are shown in Fig.^ 
with downwards pointing triangles. In the scalar 
channel the mass agrees with the measurements 
of p8j, but the vector mass is slightly smaller, 
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Figure 8. mw in the Landau gauge (after |2q]). 



likely due to the better projection to the ground 
state. In the scalar channel also standard pla- 
quette 'IF-bair operators are used together with 
the blocked S'-operators. In the symmetric phase 
very little mixing is seen between these operators; 
furthermore, the masses in the VF-ball sector are 
observed to be almost identical to the 3D SU(2) 
++ glueball masses measured at the same value 
of /? G Jp. 

In pqj mw was measured by fixing to the Lan- 
dau gauge and using Af -operators to measure 
the correlations. The results are shown in Fig.^: 
in the broken phase, the results are equal to those 
measured with the operators of type eq. (|2l|). 
However, in the symmetric phase mw ~ 0.35<7|, 
independent of m\ (and T), in strong contrast to 
the ^-operator in Fig.^. The point at k = in 
Fig.^ is a pure gauge result, i.e. calculated with- 
out the Higgs field. 

Can one understand this behaviour analyti- 
cally? In the symmetric phase the perturbative 
expansion breaks down. The system may essen- 
tially behave like 2+1D SU(2) gauge theory at 
T = 0, which is confining. The Higgs field can 
be interpreted as a scalar quark, and the phys- 
ical states are now VF-balls and (<f)'(f>) bound 
states. Using this strong coupling picture and 
the static 4>-<j) potential determined with lattice 
simulations frill , Dosch et al.j46| calculated the 
mass spectrum of the bound states analytically, 
reproducing the pattern in Fig.^ quite well. On 
the other hand, the 1-loop Schwinger-Dyson gap 
equation calculation by Buchmiiller and Philipsen 



11 



|38[ gives an approximate result mw ~ 0.28 g\ in 
the symmetric phase. This is obviously in conflict 
with the vector masses in Fig.0, but is close to 
the Landau gauge fixed result, and is shown as a 
horizontal dashed line in Fig. ||. 

The static potential and the string tension a 
in the symmetric phase have been measured in 
|16|J2S|] . Also in this case no significant depen- 
dence of ajg\ on A (m* H ) was observed. When 
m§ (T) is large, the value of a was close to the 
pure 3D SU(2) gauge theory value |iq] . At large 
distances one expects the screening behaviour to 
set in; however, this has not been observed yet in 
the distances presently allowed by the available 
resources. 

These results, while clearly supporting the con- 
finement picture in the symmetric phase, also 
show that the gauge degrees of freedom decou- 
ple almost completely from the Higgs field. In- 
deed, all the measurements which involve only 
gauge fields (plaquctte correlators, Landau gauge 
Af correlator, string tension) give results almost 
identical to the pure SU(2) gauge theory. 

5. SPHALERON TRANSITION RATE 

Due to the axial anomaly the baryon number 
is not conserved in the EW theory. In order 
to quantify the effects of the EW physics to the 
baryon number of the Universe [0 the knowledge 
of the sphaleron transition rate T is essential. In 
terms of the gauge fields, L is the diffusion rate 
of the topological charge B{t) = AAcs(t) where 
Acs is the Chern-Simons number. The effective 
potential is periodic to the Acs -direction: there 
are large gauge transformations which change 
Acs by unity. The goal is to calculate the rate 
of the dynamical processes which change Acs, 
driving the configuration from one minimum to 
a neighbouring one. Since B{t) can be described 
as a random walk in the periodic potential, 



(B 2 (t)h 



TVt, 



t 



(22) 



where (-)t is the canonical expectation value p7l . 
From analytical considerations one expects 
r = n(awT) 4 , where 



= const. 



T > T n 



fexp[-E sph (T)/T] T<T C 



(23) 
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The calculation of T is non-perturbative. Since 
r is a real-time transition rate, the standard 
imaginary time lattice formalism is not easily 
applicable. However, one expects that at high 
T the transitions occur predominantly through 
thermally activated classical processes. This sug- 
gests the following strategy: B{t) is calculated 
by solving the classical equations of motion, and 
the results are averaged over a canonical ensemble 
Q]. It is important that the Gauss constraint is 
satisfied. This is simple in the 1+1D U(l)-Higgs 
model, which has been studied in detail as a pro- 
totype model @Q, but for 3+1D SU(2) correct 
methods have been derived only recently [ plp^ ] . 

In the 1+1D U(l)-Higgs model the numerical 
and analytical results agree (however, in the high 
T phase some lattice spacing dependence remains, 
53 1). Ambj0rn and Krasnitz |54~ 



see 



measured 

T in the high T phase of the pure gauge SU(2) 
theory, with the result V = K,(awT) A , k = 1.09 ± 
0.04. The result was seen to be independent of 
(large) V and (small) a, indicating that it survives 
to the continuum limit. This is supported by the 
result by Moore |Q , who also measured the linear 
response T M of Acs to a chemical potential, r M 
'21'. 

In the 3+1D SU(2)-Higgs model it has been 
analytically estimated that K « 0.01 in the high 
T phase p5[ . Tang and Smit have recently calcu- 
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lated T in this model |3(J ; the results are shown in 
Fig.|j| In the high T phase k » 1, which agrees 
with the SU(2) result but is in strong contrast 
to the analytical estimate. At T c n decreases 
sharply, and remains 0.1-0.2 when T < T c . The 
disagreement with the analytic form in eq. (|2^) is 
dramatic: no exponential suppression is seen, and 
the rate is 10 3 -10 8 times larger than expected! 
Similar behaviour has been observed by Ambj0rn 
and Krasnitz |6). The success of the 1+1D U(l)- 
Higgs model makes this conflict even more strik- 
ing, and obviously the reason for this must be 
understood before we can have trust in the re- 
sults. 

6. CONCLUSIONS 

The numerical EW simulations have been very 
successful: to a large extent, the static thermo- 
dynamic properties of the EW phase transition 
have been "solved" . The transition is strongly 
1st order at small ran, becomes rapidly weaker 
with increasing ran and at mn ~ 80GeV (for 
SU(2)-Higgs) the line of transition ends at a crit- 
ical point, after which only a smooth cross-over 
remains. For practical purposes, the accuracy is 
good enough for most of the static quantities rele- 
vant to the transition. The properties of the criti- 
cal point itself (exact location, exponents) are not 
yet so well known, this being arguably the most 
difficult point in the phase diagram. These re- 
sults rule out the MSM baryogenesis: it is not any 
more possible that the transition is strong enough 
to produce the observed B asymmetry Q. 

The 2-loop perturbation theory yields a good 
guideline for the transition at itih < 70 GeV, al- 
though deviations are clearly seen. The 3D ef- 
fective theories and their accuracy are now fully 
understood theoretically, and the good general 
agreement with the 4D simulations is very en- 
couraging. However, more comparisons should 
be done in order to fully quantify the accuracy: 
this is especially important since the 3D theory 
provides a method for investigating realistic EW 
(+ beyond) theories without any of the problems 
usually caused by fermions, chiral or not. 

There are still unanswered questions with the 
current results, which must be addressed: the role 



of the different correlation lengths in the symmet- 
ric phase is still not fully clarified, and the con- 
tradiction between the analytical and numerical 
results for the sphaleron rate must be understood. 
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